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摘要 : 基于 经 典 极 移 理论 ， 创 新 性 地 提出 了 一 种 估计 地 球 钱 德 勒 摆动 (Chandler wobble) 周期 
和 品质 因子 的 数学 模型 和 快速 算法 ， 该 方法 对 模型 误差 的 假设 较 弱 ， 估 计 结 果 的 统计 性 质 更 优 。 
采用 1993 年 1 月 一 2009 年 12 月 时 间 段 内 的 极 移 观 测序 列 和 大 气 、 海 洋 角 动量 数据 ， 新 方法 对 
钱 德 勒 摆动 周期 和 品质 因子 的 点 估计 分 别 为 (430.8 + 0.50) d 和 62.6 士 9.63， 与 最 新 研究 结果 
近 ; 与 之 对 应 的 90% 置信 区 间 分 别 为 (430.0, 431.6) d 和 (43.5, 75.7)， 这 一 极 窗 的 置信 区 间 提 
高 了 点 估计 结果 的 可 靠 性 ， 有 助 于 更 加 精确 地 认识 地 球 钱 德 勒 摆动 的 激发 和 维持 机 理 。 
X 键 dg: 地 球 自转 变化 ， 极 移 理论 ， 角 动量 ， 优 化 算法 ， 自 助 法 
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1 引 言 


钱 德 勒 摆动 是 地 球 自转 的 一 个 本 征 模 。 自 1891 年 首次 被 发 现 以 来 ， 引 起 钱 德 勒 摆动 的 
= 物理 机 制 和 描述 它 的 两 个 关键 参数 (周期 和 品质 因子 ) 一 直 是 科学 界 关注 的 热点 问题 。 空 间 
p 大 地 测量 技术 为 检测 钱 德勤 摆动 提供 了 丰富 、 高 精度 和 高 时 空 分 辨 率 的 观测 资料 ， 人 们 利 
= 这 些 高 精度 的 观测 资料 ， 在 引起 钱 德 勒 摆动 的 物理 机 制 有 关 领 域 的 研究 取得 了 很 大 进展 ， 
例如 Smith 和 Dahlen”, Hip, Mathews A", WR Chen 等 人 " ZEB EER 
了 各 种 地 球 物理 因素 (EE HRA, WERE AE) 对 钱 德 勒 摆 动 的 贡献 ， Wahr'”、 
Haa Liao AT", DUE Fang 等 人 "研究 了 大 气 、 海 


Eubanks”, Barnes =A", Gross 
洋 等 地 表 流 体 激发 钱 德 勒 摆动 的 物理 机 制 。 尽 管 如 此 ， 由 于 理论 假设 的 局 限 性 ”"、 部 分 观测 
数据 的 不 确定 性 ””， 以 及 观测 数据 的 不 完备 性 "等 原因 ， 到 目前 为 止 ， 引 起 钱 德 勒 摆动 的 物 
时 机 制 ， 择 动 周 期 和 品质 因子 的 取 值 也 一 直 没 有 明确 的 结论 。 
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由 于 钱 德 勒 摆 动 与 很 多 地 球 物理 因素 有 关 ， 因 此 ， 对 其 精确 的 检测 具有 很 重要 的 科学 价 
值 。 目 前 ， 关 于 钱 德 勒 摆 动 周 期 和 品质 因子 的 计算 和 估计 方法 主要 分 为 如 下 四 类 : 

(1) 谱 分 析 方法 。 直 接 对 极 移 观 测序 列 进行 谱 分 析 ， 从 而 获得 钱 德勤 摆动 周期 和 品质 因 
了 由 
(2) 随机 激发 假设 方法 。 假 设 钱 德 勒 摆动 的 激发 机 制 是 随机 的 ， 采 用 统计 方法 来 估计 钱 
德勤 摆动 周期 和 品质 因子 ” ”。 

(3) 地 球 物理 激发 理论 方法 。 采 用 实际 观测 数据 ， 研 究 大 气 、 海 洋 以 及 陆地 水 等 地 球 物 
理 机 制 对 钱 德勤 摆动 的 影响 "” 79 

(4) 半 解 析 方法 。 基 于 一 定 的 理论 模型 和 实际 观测 资料 ， 例 如 海洋 、 地 慢 、 核 慢 耦 合 和 
三 轴 椭 球 等 等 ， 给 出 钱 德 勒 摆 动 周 斯 和 品质 因子 的 解析 表达 式 ”“。 
本 文 将 前 三 类 方法 获得 的 周期 和 品质 因子 称 之 为 观测 值 ， 将 最 后 一 类 方法 获得 的 称 之 
为 理论 值 。 表 1 统计 了 现 有 方法 对 钱 德 勒 摆 动 周 期 和 品质 因子 的 计算 和 估计 结果 。 从 表 1 可 
以 看 出 ， 品 质 因 子 的 参数 估 值 具有 极 大 的 取 值 范围 (36, 1000)， 尽 管 最 新 研究 已 经 将 其 置信 
区 间 缩 小 至 (56, 255)"”， 但 仍然 比较 大 。 


表 1 钱 德 勒 摆动 周期 和 品质 因子 的 点 估计 和 区 间 估 计 的 现 有 结论 ” 


方法 周期 /太阳 品质 因子 信和 度 ” 时 间 跨 度 /a 参考 文献 
429.73 土 7.09 一 一 80 17 
谱 分 析 法 432.71 士 0.42 36.7 E 61 18 
4340.5 = = 106 19 
433.2:-3.4 61 (37, 193) lo 68 20 
434.8+2.0 96 (50, 300) lo 76 21 
随机 激发 假设 433.343.1 170 (47, 1 000) lo 78 22 
433.041.1 179 (74, 789) lo 86 23 
433.141.7 一 一 147 24 
439.5 士 2.1 72 (30, 500) 90% 9 26 
433.7£1.8 49 (35, 100) lo 11 25 
地 球 物理 激发 理论 429.4 107 = 10 27 
432.98 97 = 60 28 
f 430.9 土 0.7 127 (56, 255) 95% 25 16 
i 430.84+0.50 62.6 (43.5, 75.7) 90% 17 本 工作 
426.73-2.6 — — — [1] 
434.8 71 一 一 [3, 4] 
半 解 析 法 430.3 88.4 — 20 [5] 
433.03 100.20 一 20 [29] 


[13 


注 : @ 该 表格 参考 了 Gross 的 表 12 和 Nastula 等 人 W 的 表 2， 并 引用 了 更 多 文献 ，@ 这 一 列表 示 品 质 
因子 的 置信 水 平 ， 其 中 1o 表示 单位 标准 差 置信 区 间 荆 2 。 


在 前 人 工作 基础 上 ， 本 文 创新 性 地 提出 了 一 种 用 于 估计 钱 德 勒 摆动 周期 和 品质 因子 的 
数学 模型 和 快速 算法 。 新 方法 对 模型 误差 的 假设 更 弱 ， 估 计 结 果 的 统计 性 质 更 优 ， 新 方法 还 
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采用 了 机 制 明确 的 优化 算法 ， 计 算 效 率 更 高 。 


2 极 移 理论 


假设 地 球 不 受 外 部 天 体 的 引力 作用 ， 在 地 固 坐 标 系 下 ， 角 动量 定理 可 以 写成 刘 维 尔 方 
程 ， 即 : 


L 
C twexb=0, (1) 


其 中 ，w 是 地 固 系 相 对 惯性 系 的 转动 角速度 。 用 1 表示 地 球 自转 的 平均 角速度 ， 则 w 可 以 


w = (Q(m;,m»,1- ma) . (2) 
CL 是 地 球 自转 角 动量 ， 在 地 固 系 下 ， 其 表达 式 是 


bed gh. (3) 


其 中 ，h = (h, ho, hg)? 是 地 球 的 相对 角 动 量 ， 了 是 地 球 的 惯性 张 量 ， 它 的 表达 式 为 


4 0 0 Alı Ah» Ahs 

[> 

I = 0 4 0 + Al Aloo Als , (4) 
00 C 


Kp, AFC 是 地 球 的 主 惯性 矩 。 将 diag( A, A, C) 作为 T 的 初始 值 或 无 摄 动 项 ， 则 上 式 
第 二 部 分 代表 的 就 是 了 的 摄 动 ， 它 是 随时 间 变 化 的 。 

根据 Wahr” 的 公式 (2.17), (2.21) 和 (2.24)， 可 知 极 移 m(t) 与 激发 函数 L(t) 满足 如 下 
微分 方程 : 


zit) + mt) =I SAE ， (5 


其 中 ，A EIRE, A 与 Am 的 关系 为 4 = 1.128434, BRA = 1.1294, 5 Cow 
是 钱 德 勒 频率 ， 在 这 里 它 只 是 一 个 实数 ; mt) 的 表达 式 为 m(t) = ma (t) + ima(t); v(t) 的 
AX y(t) = rac ay ne )-il(t)), Wt) 是 包含 了 大 气 、 海 洋 等 因素 的 角 
动量 ， 它 的 表达 式 为 : 


H(t) = h(t) + QAT(t) ， (6) 


RP, A(t) = hi(t) 十 jh2(t)， 即 相对 角 动 量 ， 由 于 它 与 物质 相对 于 地 固 系 的 运动 速度 有 关 ， 
所 以 又 称 之 为 运动 项 ;AT(t) = ATs(t) 十 iA1z3(t)， 即 地 球 惯性 张 量 的 摄 动 ， 由 于 它 只 与 物 
质 的 密度 分 布 有 关 ， 所 以 又 将 QAT(t) 称 之 为 质量 项 。 

考虑 地 慢 的 滞 弹 性 和 海洋 的 非 平衡 响应 ， 可 以 先 将 式 (5) 中 的 实数 oc 车 换 为 复数 


pu (+ 如 ) 其 中 到 和 Q. 分 别 代表 钱 德 勒 摆动 周期 和 品质 因子 ， 再 对 其 化 简 ， 结 


Te 2Q. 


-H 
H 


年 期 
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RA 
: i i. 
m(t) —io.m(t) = m (mo = ait) i (7) 


这 里 对 式 (7) 做 一 点 说 明 ， 传 统 处 理 方法 是 先 对 式 (5) 化 简 ， 再 将 其 中 的 oow PION oe, 
即 先 化 简 再 蔡 换 ， 如 文献 [31] 中 的 式 (13) 和 (14)， 而 本 工作 则 是 先 普 换 再 化 简 。 由 于 这 两 
种 处 理 方法 对 m(t) 的 影响 不 超过 1%"， 并 且 观 测 数据 的 测量 精度 有 限 ， 所 以 可 以 认为 这 两 
种 做 法 一 致 。 考 虑 极 移 m(t) 与 观测 极 移 p(t) 的 关系 ， 可 得 ”: 


m(t) = ptt) - Zale) ， (8) 


其 中 ，p(t) = a(t) — iy(t), x(t) 和 y(t) 就 是 由 国际 地 球 自转 服务 (International Earth 
Rotation and Reference Systems Service, IERS) 或 其 他 组 织 和 机 构 提供 的 极 移 观 测序 列 ， 其 
中 y(t) 以 指向 90° W 为 正 。 将 式 (8) 代入 到 式 (7) 中 ， 可 得 


Blt) = ioli) = ——1 (9) 


对 这 个 微分 方程 进行 求解 ， 其 中 p(t) 的 自由 项 在 经 过 很 长 一 段 时 间 后 逐渐 衰减 消失 ， 只 剩 
下 受 迫 激发 项 ， 因 此 p(t) 的 解析 表达 式 为 : 


t 
p(t) = -ent f eice77T(r)dr . (10) 


3 估计 算法 


基于 式 (10)， 本 章 对 误差 进行 基本 假设 ， 提 出 用 于 估计 T. AQ. 的 数学 模型 ， 同 时 ， 基 
于 这 些 误差 假设 和 新 的 数学 模型 ， 引 入 T. 和 Q 的 区 间 估 计 方 法 一 一 自助 法 。 
3.1. fait T, 和 Q. 的 快速 算法 

假设 数据 是 等 间隔 采样 的 ， 采 样 间隔 为 比 ， 数 据 长 度 为 W， 初 始 历 元 为 ta。 根据 式 
(10) 可 以 给 出 离散 数据 满足 如 下 一 阶 自 回归 形式 : 


: tn41 
, i . : 
eie hp(t, 41) u p(ts) ne E" eicetn f e i7e7 II (7)dr 
m tn 


(t, =to +nét , n=0,1,2---,N—1) . (11) 
考虑 到 II(t) 不 完全 是 大 气 和 海洋 的 角 动 量 ， 可 能 还 包含 了 不 可 观测 的 和 未 知 的 因素 ， 如 在 


式 (6) AY A(t) RI AT(t) 中 ， 它 们 不 仅 包 含 了 地 表 流 体 的 作用 ， 还 包括 了 地 球 内 部 的 作用 (如 
地 磁 急 变 ”)， 而 后 者 不 可 直接 观测 ”"。 所 以 ， 可 以 将 I(t) 分 解 为 如 下 三 个 部 分 : 


H(t) = Hovy(t) + Hos (t) + Tros(t) ， (12) 


= Hows (t) + Her (t) IWR A, Hon (0) 是 与 之 对 应 的 实际 观测 数据 ，Zer. 人 
观测 数据 的 测量 误差 (BMR); Mres (t) 是 上 面 提 到 的 其 他 不 可 观测 的 或 未 知 的 激发 
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源 。 这 里 将 式 (12) 的 后 两 项 合并 ， 记 为 [oor(t)， 简 称 模型 误差 。 最 后 ， 假 设 模型 误差 对 应 
的 离散 序列 Tus (5) 是 一 个 均值 函数 为 0、 一 阶 自 相 关 函 数 为 0 的 随机 噪声 。 

Hi E) = Holt) + Hass (0) RAR (11) 中 ， 对 关于 Molt) 的 积分 用 辛普森 公式 展开 ， 
THAF nalt) 的 积分 用 梯形 公式 展开 ， 并 引入 以 下 表达 式 : 


bı = exp(—idto,/2) 


b2 = exp(—idtoc) (13) 
EN 


Xn, = —— How (tn+0.5) 
X2 = Ppltn+1) + S4 Ho (tnt) : (14) 
Ya = p(tn) I S4 How (ts) 
利用 式 (13) 和 (14) 将 式 (11) 展开 并 化 简 为 : 
Yn = Xn Pi = X522 = En + En+1b2 , (15) 


其 中 , n=0,1,2,---,N—-2. SL AAT Y, Xi, Xo, e Me, X (15) 可 以 被 改写 为 回归 方程 
的 形式 : 


Y — X16, — X505 = E + €f» é (16) 
考虑 到 序列 Mme lta) 的 一 阶 自 相关 函数 为 零 ， 因 此 eNé = 0 成 立 ( 上 标 瑟 代表 矢量 的 共 

斩 转 置 )， 可 以 获得 残 差 平方 和 的 解析 表达 式 ; 
IY — X181 — X222? 


llel? = 1+ BEB, (17) 
其 中 ，|| .| 代表 矢量 的 三 范 数 。 为 了 求解 61， 可 以 将 式 (17) 转化 为 约束 优化 问题 ， 即 : 
nllell , s.t.6 = 67 . (18) 
使 用 拉 格 朗 日 乘 子 法 求解 ， 定 义 拉 格 朗 日 函数 L(B1, 82, A) 为 : 
L(81, B2, A) = |lel — A(B2 — Br) . (19) 
Xf L(B1, 82, A) 求 3 个 偏 导数 ， 令 它们 都 为 0， 再 消去 入 可 得 : 
(Y - X. - iX)" [22 (X2 + BY) + (1 |B )X)=0 . (20) 


将 等 式 (20) 的 左边 记 为 F(p RERE f(8) 的 零点 即 可 。 采 用 数值 分 析 中 的 双 点 制 线 法 过 
RRA, ARRU F: 
JU AB) ( 1 


2 Re(8* — 8*-1)  Im(f* 一 =) TUE? ， (21) 


ghe = m 


C— 
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其 中 ，Re() 和 Im() 分 别 代表 取 虚 数 的 实 部 和 虚 部 。 此 时 ， 只 需 给 定 初 什 89 = exp [77 
(itag) | Bae be 其中 < 是 一 个 比较 小 的 复 随机 数 ， 将 四 和 中 代入 式 (21)， 反 


Ak. WR BP 与 BF 的 差距 小 于 某 一 阐 值 (这 里 取 为 10-4)， 则 跳出 循环 并 认为 
B — Br, WIE Ê DWT ARH Ê, 和 0。。 

总 之 ， 该 方法 采用 了 一 种 机 制 明确 的 算法 来 直接 获得 最 优 的 T. 和 @。， 而 不 用 在 某 个 特 
定 区 域内 比较 每 个 网 格 点 上 的 结果 ， 例 如 Nastula 和 Gross “采取 的 方法 (下 面 将 其 简称 为 
网 格 搜 索 法 )， 因 此 本 方法 具有 更 高 的 计算 效率 。 
3.2 ”区 间 估 计 的 方法 一 一 自助 法 
昌 前 关于 钱 德 勒 摆动 周期 和 品质 因子 的 区 间 估 计 都 是 基于 蒙特 卡 洛 方法 开展 的 ””， 这 
种 方法 需要 事先 获取 模型 误差 melt) 的 理论 分 布 。 根 据 式 (12) 可 知 : = Tres (t) 是 未 知 的 或 
不 可 观测 的 ， 假 设 它 服从 任何 一 种 分 布 可 能 都 不 合理 ， 此 外 ， 随 着 观测 水 平 的 提高 ， 观 测 数 
据 的 测量 精度 也 在 逐渐 提高 ， 因 此 观测 误差 Melt) 具有 异 方差 性 。 综 合 这 两 个 因素 可 以 发 
现 ， 理 论 模 型 难以 描述 Tae(h 的 分 布 函数 ， 因 此 蒙特 卡 洛 方法 的 区 间 估 计 结 果 可 能 不 够 精 
傅 。 


基于 此 ， 本 工作 采用 自助 法 (Bootstrapping) 做 区 间 估 计 ， 该 方法 的 优势 是 不 需要 知道 
IL (t) 的 分 布 函 数 。 自 助 法 是 一 种 现代 的 非 参 数 统计 方法 ， 最 早 由 Efron” 提 出 ， 经 过 近 
几 十 年 的 发 展 ， 它 已 经 具备 了 扎实 的 理论 背景 并 拥有 广泛 的 应 用 范围 ”。 自 助 法 的 操作 流 
程 简捷 ， 其 体 为 : (1) 将 包含 到 个 样品 的 初始 数据 集 当做 总 体 ， 每 次 从 中 有 放 回 地 抽取 C 
个 样品 ， 如 此 重复 B 次 ， 就 得 到 了 B 个 新 样本 ; (2) 对 每 个 再 抽样 的 样本 分 别 做 参数 估计 ， 
可 以 得 到 每 个 参数 的 B 个 估 值 ，(3) 基于 这 B 组 估 值 ， 可 以 给 出 基于 初始 数据 集 的 区 间 估 
计 ， 甚 至 研究 参数 估 值 的 有 效 性 。 

本 工作 采用 自助 法 来 对 工 和 Q. 进行 区 间 估 计 ， 按 照 以 上 步骤 将 式 (15) 中 的 W 一 工 个 
方程 当成 总 体 (这 一 做 法 在 参考 文献 33] 中 有 详细 的 描述 )， 每 次 从 中 有 放 回 地 抽取 N 一 1 
个 样品 ， 然 后 根据 式 (16) 一 (21) 对 再 抽样 的 样本 做 参数 估计 ， 如 此 重复 多 次 ， 可 以 获取 多 组 
T. MQ. 的 参数 估 值 ， 据 此 给 出 到 Al Qe 的 区 间 估 计 结 果 。 


O 4 数据 及 预 处 理 


对 于 极 移 观测 序列 ， 我 们 采用 Ratcliff 和 Gross” 解 算 的 COMB2018 数据 集 ， 该 数据 
集 在 1990 年 后 的 时 间 段 具有 非常 高 的 精度 。 大 气 角 动量 (atmospheric angular momentum, 
AAM) 数据 来 自 美 国 国 家 环境 预报 中 心 / 大 气 研究 中 心 (National Centers for Environmental 
Prediction / National Center for Atmospheric Research, NCEP / NCAR) 数据 集 ”， 该 数 
据 集 考虑 了 地 形 因素 的 影响 ， 精 度 更 高 ""。 海 洋 角 动量 (oceanic angular momentum, OAM) 
数据 来 自 TERS 特殊 海洋 局 (Special Bureau for Ocean) 的 ECCO-kf080i 数据 集 。 这 三 个 数 
据 集 的 采样 间隔 都 是 1d， 选 取 的 公共 时 间 跨 度 为 1993/01 一 2018/12。 在 参数 估计 之 前 ， 对 
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这 些 数 据 做 如 下 预 处 理 ”“” : 

(1) 扣除 三 个 序列 中 的 周年 项 、 季 节 项 和 线性 趋势 项 。 采 用 最 小 二 乘法 从 原始 序列 中 扣 
除 频率 为 1, 2, 3, 4, 5, 6, 7, 8 cpy 的 谐 波 信号 ， 两 个 周期 分 别 为 13.661 d 和 13.633 d 的 潮汐 
项 ， 一 个 一 阶 多 项 式 函 数 。 

(2) 扣除 低频 项 。 除 线性 趋势 项 外 ， 极 移 序列 中 的 低频 信号 一 般 被 认为 来 自 地 球 内 部 ， 
而 这 些 激发 源 不 能 被 直接 观测 ， 据 此 设计 一 个 高 通 滤波 器 ， 过 滤 掉 三 个 序列 中 周期 大 于 2 a 
的 频段 。 

预 处 理 后 的 数据 展示 在 图 1 中 ， 从 该 图 可 以 看 出 : 钱 德 勒 摆动 的 振幅 一 直 处 于 调制 中 ， 
总 体 在 减 小 ， 而 在 2010 一 2019 年 期 间 ， 钱 德勤 摆动 的 振幅 已 经 减 小 到 非常 显著 的 程度 (这 
也 可 以 在 Wang 等 人 "的 图 8 中 得 到 验证 )， 对 于 预 处 理 后 的 角 动量 序列 ， 与 平稳 随机 噪声 
非常 相似 ， 并 且 实 部 振幅 比 虚 部 振幅 稍 小 。 


钱 德 勒 摆动 的 实 部 l 钱 德 勒 摆动 的 虚 部 
200 4 | 
g 100 4 100} 
B g 
pg 0 只 0 
b x E» 
1 
S -100H 3$ -100 
-200r ] -200 
1995 2000 2005 2010 2015 1995 2000 2005 2010 2015 
时 间 / 年 份 时 间 / 年 份 
x1024 AAM+OAM 的 实 部 l xl024 ”AAM-+OAM 的 虚 部 
n 
a 
Z 
= 
pE 
E 


1995 2000 2005 2010 2015 1995 2000 2005 2010 2015 
时 间 / 年 份 时 间 / 年 份 


TE: 第 一 行 是 预 处 理 后 的 极 移 序 列 p 中 的 Chandler 项 ， 第 二 行 是 预 处 理 后 的 角 动 量 序列 Tobv， 是 大 气 角 动 
(AAM) 和 海洋 角 动 量 (OAM) 加 和 后 的 结果 。 


1 ， 预 处 理 后 的 数据 序列 图 


5 结 d 


为 了 更 全 面 地 评估 新 方法 的 稳定 性 和 可 靠 性 ， 本 工作 选用 三 个 不 同 的 时 间 段 对 研究 结 
果 进 行 分 析 和 讨论 ， 分 别 是 1993 一 2010 年 、2000 一 2019 年 和 1993—2019 年 。 在 每 个 时 间 
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段 都 使 用 自助 法 做 1000 次 再 抽样 ， 每 个 时 间 段 都 得 到 1000 28 T, A Qe. 3: 3000 组 估计 值 
的 统计 直方 图 如 图 2 所 示 ， 它 们 的 均值 、 方 差 和 9096 置信 区 间 如 表 2 所 示 。 


150 200 
1993 一 2010 年 1993 一 2010 年 


£ 
= 100 
Acs 
50 
433 434 0 80 100 120 140 
2000—20194F 2000—20194F 
100 
< 
3 
= 50 
fy 430 431 432 433 484 80 100 120 140 


1993—20194 


1993—20194F 


计数 /个 


429 430 431 432 433 40 60 80 100 120 140 
周期 /d 衰减 因子 


2 3000 组 周期 和 品质 因子 的 参数 估计 值 的 统计 直方 图 


表 2 3000 组 周期 和 品质 因子 的 参数 估计 的 均值 、 方 差 和 90% 置信 区 间 
HE 周期 /d 品质 因子 
点 估计 ”标准 差 90% 置信 区 间 | 点 估计 标准 差 90% 置 信 区 间 
430.8 0.50 (430.0, 431.6) | 62.6 9.63 (43.5, 75.7) 
432.2 0.68 (431.1, 433.3) | 76.2 25.4 (24.3, 97.7) 
430.7 0.46 (429.9, 431.5) | 52.5 6.56 (40.3 60.8) 


1993 一 2010 4 
2000 一 2019 4 
1993—2019 4 


[un [un [un 


从 图 2 和 表 2 的 统计 结果 可 以 看 出 ， 第 一 个 时 间 段 (1993 — 2010 年 ) 和 第 三 个 时 间 段 
(1993 一 2019 4E) 给 出 的 到 AQ. 的 区 间 估 计 几 乎 是 重合 的 ， 并 且 A 的 方差 很 小 ， 对 于 Q 
的 估计 ， 其 对 数据 的 选择 比较 敏感 ， 这 也 是 现 有 估计 方法 共同 的 问题 (如 表 1 所 示 )。 对 比 表 
1 和 表 2 可 以 看 出 ， 新 方法 对 Q. 的 估计 具有 更 小 的 方差 ， 显 示 了 更 优 的 统计 有 效 性 。 
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从 图 2 和 表 2 的 统计 结果 还 可 以 发 现 ， 第 二 个 时 间 段 (2000 一 2019 年 ) 与 其 他 两 个 时 间 
段 的 参数 估计 结果 差距 较 大 。 这 可 能 与 近年 来 钱 德勤 摆动 的 振幅 逐渐 减 小 有 关 ， 该 时 间 段 的 
估计 结果 可 能 具有 一 定 的 系统 性 偏差 ， 我 们 将 在 未 来 的 研究 中 进一步 探讨 。 

综合 以 上 结果 ， 可 以 确定 第 一 个 时 间 段 (1993—2010 4E) 的 参数 估计 更 加 合理 ， 即 区 
和 Qe 的 点 估计 为 430.8 d 和 62.6， 两 者 的 90% 置信 区 间 分 别 为 (430.0, 431.6) d 和 (43.5, 
75.7)。 在 该 时 间 段 内 的 参数 置信 区 间 更 宽 ， 具 有 更 大 的 概率 包含 Qe 的 真 值 ， 是 一 个 比较 保 
守 的 估计 。 此 外 ， 这 一 结果 与 Mathews 等 人 "及 Nastula 和 Gross 的 结果 很 接近 ， 再 一 次 
验证 了 新 方法 的 可 靠 性 。 


6 总 结 与 讨论 


本 文 从 经 典 的 极 移 理论 出 发 ， 给 出 了 一 种 用 于 估计 钱 德 勒 摆动 周期 和 品质 因子 的 新 方 
法 。 新 方法 对 于 误差 的 假设 较 弱 ， 并 采用 更 加 合理 的 自助 法 进行 区 间 估 计 ， 参 数 估 计 的 结果 
\ 有 更 优 的 统计 性 质 。 此 外 ， 新 的 点 估计 算法 采用 了 机 制 明 确 的 双 点 割 线 法 来 直接 获取 最 优 
解 ， 与 传统 的 网 格 搜索 法 相 比 具有 更 高 的 计算 效率 。 

对 地 球 钱 德勤 摆动 周期 和 品质 因子 的 最 优点 估计 结果 分 别 为 : 430.8 d 和 62.6， 两 者 的 
标准 差 为 0.50 d 和 9.63， 最 优 的 区 间 估 计 分 别 为 (430.0, 431.6) d 和 (43.5, 75.7)。 这 一 结果 
与 Mathews 等 人 "及 Nastula 和 Gross’ 的 结果 很 接近 ， 且 我 们 的 参数 估 值 具有 更 小 的 方 
差 ， 再 次 验证 了 新 方法 的 可 靠 性 和 优越 性 。 
由 于 目前 只 采用 了 大 气 和 海洋 的 角 动 量 数据 ， 对 地 球 钱 德 勒 摆动 周期 和 品质 因子 的 最 
优 估计 可 能 还 存在 一 定 偏差 。 未 来 我 们 将 考虑 陆地 水 或 全 球 水 文 角 动 量 的 影响 ， 进 一 步 提高 
钱 德 勒 摆动 周期 和 品质 因子 的 估计 精度 。 此 外 ， 随 着 对 钱 德 勒 摆动 认识 的 逐步 提高 以 及 观测 
数据 的 精度 提升 ， 本 工作 提出 的 新 方法 将 对 地 球 钱 德 勒 摆动 的 精准 预报 提供 可 能 。 
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Abstract: Based on the classic polar-motion theory, we establish a new method contains 
mathematical model and fast algorithm, to estimate the period and Q Value of the Earth's 
Chandler wobble. This method makes weak assumptions on the model error, and the esti- 
mated results have better statistical properties. Given the polar-motion, atmospheric and o- 
ceanic angular momentum observations during 1993/01 —2009/12, we apply the new method 
to estimate the Chandler wobble's period and Q value, the preferred point estimations are 
430.8( 士 0.50) solar days and 62.64( 士 9.63) in respective, which are close to the latest re- 


search results. Moreover, the 90% confidence interval results are (430.0, 431.6) solar days 


and (43.50, 75.66) respectively, such a extremely narrow confidence interval can increase the 
reliability of the point estimations, which is helpful for understanding of the excitation and 


maintenance mechanism of the Chandler wobble more precisely. 


Key words: Earth's rotation variation; polar-motion theory; angular momentum; optimiza- 


tion algorithm; bootstrapping 


